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Abstract 

Progressive failure analyses results are pre- 
sented for composite panels with and without a 
cutout and are subjected to in-plane shear load- 
ing and compression loading well into their post- 
buckling regime. Ply damage modes such as ma- 
trix cracking, fiber-matrix shear, and fiber fail- 
ure are modeled by degrading the material prop- 
erties. Results from finite element analyses are 
compared with experimental data. Good agree- 
ment between experimental data and numerical 
results are observed for most structural configu- 
rations when initial geometric imperfections are 
appropriately modeled. 
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The use of composite materials for aircraft 
primary structures can result in significant ben- 
efits on aircraft structural cost and performance. 
Such applications of composites materials are ex- 
pected to result in a 30-40 percent weight savings 
and a 10-30 percent cost reduction compared to 
conventional metallic structures. However, un- 
like conventional metallic materials, composite 
structures fail under different failure modes such 
as matrix cracking, fiber-matrix shear failure, 
fiber failure, and delamination. The initiation 
of damage in a composite laminate occurs when 
a single ply or part of the ply in the laminate fails 
in any of these failure modes over a certain area 
of the structure. The initiation of damage does 
not mean that the structure cannot carry any 
additional load. The residual load bearing capa- 
bility of the composite structure from the onset 
of material failure or initiation of damage to fi- 
nal failure can be quite significant. It is at the 
final failure load that the structure cannot carry 
any further load. This may be due to the fact 
that some failure modes may be benign not to 
degrade the performance of the overall structure 
significantly. Accurate determination of failure 
modes and their progression helps either to de- 
vise structural features for damage containment 
or to define fail-safe criteria. Therefore it is im- 
portant to understand the damage progression in 
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composite structures subjected to different load- 
ing conditions. 

Considerable work has been performed on 
this subject. In 1987, Talreja ([1]), Allen et al. 
([2]), and Chang and Chang ([3]), independently 
proposed progressive failure models that describe 
the accumulation of damage in a composite lam- 
inate by a field of internal state variables. The 
damage model proposed by Chang and Chang 
([3]) for notched laminate loaded in tension ac- 
counts for all of the possible failure modes ex- 
cept delamination. Chang and Lessard ([4]) later 
investigated the damage tolerance of composite 
laminates subjected to compression. Davila et 
al. ([5]) extended Chang and Lessard method 
from two dimensional membrane effects to shell- 
based analysis that includes bending. Shell-based 
progressive failure analyses that apply a mate- 
rial degradation model at every material point 
in every ply in the laminate have been reported 
in References [6] through [9]. Recently a shell- 
based progressive failure analysis which consid- 
ered large rotations based on a total Lagrangian 
approach method was presented in Reference 
[10]. 

The objective of the present paper is to de- 
velop and validate an efficient methodology that 
can predict the ultimate strength of compos- 
ite panels by taking into account ply damage 
modes and geometrical non-linear response. Re- 
sults from progressive failure analyses of compos- 
ite panels with and without a cutout subjected to 
shear and compressive loads are compared with 
experimental results. In this paper progressive 
failure results are compared with experimental 
results for fiat panels with and without cutout 
subjected in-plane shear loading. Progressive 
failure results are also presented for curved panels 
with and without a cutout and subjected to ax- 
ial compression. Although results from such ex- 
perimental studies have been compared by other 
authors with postbuckling analyses results, pro- 
gressive failure analyses for nonlinearly deformed 
structures is not reported in the open literature. 
Thus another objective for this paper is to pro- 
vide such a comparison. 

Failure Analysis 

Failure modes in laminated composite panels 
are strongly dependent on ply orientation, load- 
ing direction and panel geometry. There are four 
basic modes of failure that occur in a laminated 
composite structure. These failure modes are; 
matrix cracking, fiber-matrix shear failure, fiber 


failure, and delamination. Delamination failure, 
however, is not included in the present studies. In 
order to simulate damage growth accurately, the 
failure analysis must be able to predict the fail- 
ure mode in each ply and apply the corresponding 
reduction in material stiffnesses. The failure cri- 
teria included in the present analyses are those 
proposed by Hashin [11] and are summarized be- 
low T . 

• Matrix failure in tension and compression 
occurs due to a combination of transverse, 
<j 2 2 i an d shear stress, T\ 2 , r 13 and 723 . The 
failure index can be defined in terms of 
these stresses and the strength parame- 
ters Y and shear allowables S c . Failure 
occurs when the index exceeds unity. As- 
suming linear elastic response, the failure 
index has the form: 



where Y t is the strength perpendicular to 
the fiber direction in tension, Y c is the 
strength perpendicular to the fiber direc- 
tion in compression, and 5 c i2> 5ci3> an d 
S C 23 are the in-plane shear, and transverse 
shear strengths, respectively. 

• Fiber-matrix shear failure occurs due to a 
combination of axial stress (an) and the 
shear stresses. The failure criterion has 
the form: 



+ f° r > 0 
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where X t is the strength along the fiber 
direction in tension, and X c is the strength 
along the fiber direction in compression. 

• Fiber failure occurs due to tension or com- 
pression independent of the other stress 
component. In compression the fiber fails 
by buckling. The failure criterion has the 
form: 

e/ = -rr- for < 0 (5) 

X c 

and 

e/ = for a n > o (6) 

Xt 

To simulate the above failure modes, the 
elastic properties are made to be linearly depen- 
dent on three field variables, FVl through FV 3. 
The first field variable represents the matrix fail- 
ure, the second the fiber-matrix shearing failure, 
and the third the fiber buckling failure. The val- 
ues of the field variables are set to zero in the 
undamaged state. After a failure index has ex- 
ceeded 1.0, the associated user-defined field vari- 
able are set to 1.0. The associated field vari- 
able then continues to have the value of 1.0, even 
though the stresses may reduce to values lower 
than the failure stresses of the material. This pro- 
cedure ensures that the damaged material does 
not heal. The mechanical properties in the dam- 
aged area are reduced appropriately, according to 
the property degradation model defined in Table 
1. For example, when the matrix failure crite- 
rion takes the value of 1.0, then by the interpola- 
tion rule defined in Table 1, the transverse shear 
modulus (E y ) and the Poisson ratio (i/ 12 ) are set 
equal to zero. The field variables can be made to 
transit from 0 (undamaged) to 1 (fully damaged) 
instantaneously. Chang and Lessard’s degrada- 
tion model [4] is used in the present study. 

The finite element implementation of this 
progressive failure analysis was developed for 
the ABAQUS structural analysis program using 
the USDFLD user-written subroutine [12, 13]. 
ABAQUS calls this USDFLD subroutine at all 
material points of elements that have material 
properties defined in terms of the field variables. 
The subroutine provides access points to a num- 
ber of variables such as stresses, strains, mate- 
rial orientation, current load step, and material 
name, all of which can be used to compute the 
field variables. Stresses and strains are computed 
at each incremental load step and evaluated by 


the failure criteria to determine the occurrence of 
failure and the mode of failure. 

Numerical Examples 

To assess the predictive capability of the 
present failure analysis method, several panels 
have been analyzed and these results were com- 
pared with experimental results. Results are pre- 
sented for unstiffened panels with and without 
cutouts and a stiffened panel. The unstiffened 
panel cases are a flat panel loaded in shear and 
a curved panel loaded in compression. Results 
are also presented for a bead-stiffened panel sub- 
jected to in-plane shear loading. 

Flat Panel Loaded in In-plane Shear 

The panels subjected to pure in-plane shear 
loading were loaded using a picture frame fix- 
ture. The test sections of the bead-stiffened and 
unstiffened panel were 12.0 by 12.0 in. in size 
and the members of the picture frame were 2.75- 
in. wide and 6.75-in thick. The fixture is made 
of steel. Figure 1 shows a schematic diagram of 
a picture frame test fixture. In the finite ele- 
ment model, nodes on each member were con- 
strained for the out-of-plane displacement. Pin 
joint consists of two co-incident nodes tied in 
a multi-point constraint at the four corners of 
the panel. The displacements of the dependent 
node is made the same as that of the indepen- 
dent node, but the rotations of the co-incident 
nodes are excluded in the multi-point constraints. 
The independent node diagonally opposite to the 
loading pin is constrained for axial and transverse 
displacements. At the loading pin, applied dis- 
placement equal in magnitude in the axial and 
transverse directions at the independent node 
simulates the loading condition. The test section 
is modeled using ABAQUS four node, reduced 
integration, shear deformable S4R element [13]. 
The members of the picture frame are modeled 
using ABAQUS four node shear deformable S4 
element [13]. 

The flat unstiffened panel has a lami- 
nate stacking sequence of [±45/0/90]2$, with 
a ply thickness of 0.0056-in. and is made 
of graphite epoxy. The mechanical properties 
for the material are En — 18.5 Msi, £*>2 = 1*64 
Msi, Gi2=Gi3=0.87 Msi, G-23 =z 0*65 Msi, and 
i/ l2 = 0.3. The strength allowables are A^=232.75 
ksi, X c =210.0 ksi, If = 14.7 ksi, i;=28.7 ksi, 
S c i2=29.75 Ksi, and 5 c i 3 =5 c2 3=4.8 ksi. Experi- 
mental results for this test panel are reported in 
Reference [14], 
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A finite element model of the panel is shown 
in Figure 2. This model consists of 3425 nodes 
and 3300 elements. An imperfection based on 
static analysis results for a pressure load was 
added to the model to simulate an imperfect 
shape similar to that of a bubble (one half wave 
in each direction of the panel). Progressive fail- 
ure analysis (PFA) was carried out for this case 
with a maximum imperfection magnitude equal 
to 5% of the laminate thickness. Three integra- 
tion points through each ply thickness are used 
in the analysis for computation of section prop- 
erties. A post-buckling analysis of the panel 
with the same level of imperfection was also per- 
formed. 

The results for the flat panel loaded in shear 
are shown in Figure 3, where the load is plotted 
versus the strain normal to the fiber direction 
(£ 22 ) in the top and bottom ply (45° ply) at the 
center of the test specimen. The dashed lines 
marked FV1, FV2 and FV3 indicate the load 
level at which damage described by field vari- 
able FV 1 through FV 3 are initiated well into 
the post buckling regime. These failures are due 
to severe bending in the region. The thick solid 
line represents the experimental results and the 
thin solid line represents the postbuckling analy- 
sis results. The solid and open triangles are an- 
alytical results for the panel response that indi- 
cates progression of failure. The analysis results 
are in good agreement with the experimental re- 
sults. The final failure load obtained from the 
experiment is 54.81 kips, which is 6 % more than 
the final failure load of approximately 51.5 kips 
obtained from progressive failure analysis. Dam- 
age initiation starts as matrix cracking (FVT) at 
a load level of 37.76 kips. Fiber-matrix shear 
(FV 2), and fiber (FV 3) failure are initiated at 
the same load level of 45.37 kips. This indicates 
that the structure can carry an additional 17.04 
kips (about 30% of the final experimental load) 
after matrix cracking has been initiated. Even af- 
ter the initiation of fiber-matrix shear and fiber 
failure, the panel continues to carry an additional 
load of about 9.4 kips ( 17% of the final experi- 
mental load). The post-buckling analysis results 
diverge from the progressive failure results at a 
load level that is slightly higher than the load at 
which initiation of fiber-matrix shear and fiber 
failure occurs. At that load level the analysis pre- 
dicts a significant amount of damage that could 
have led to a considerable loss of panel stiffness. 

All the failure modes initiated near the steel 
supporting fixture close to the region along the 


diagonal which is in the panel loading direction 
and propagate in the region close to the diagonal 
and along the loading direction. Figure 4 shows 
a fringe plot of matrix cracking (FV 1) in the top 
ply after the final failure load. The dark contours 
denote failed regions of the panel. Other damage 
modes accumulate in approximately the same re- 
gion as depicted by the plot. These damage loca- 
tions are consistent with observations from exper- 
iments and do not involve delamination failure. 

Curved Panel Loaded in Axial Compression 

The curved panel has a laminate stacking 
sequence of [±45/0/90] 3 s, with a ply thickness 
of 0.005-in. The mechanical properties of the 
material used for the panel are Fn = 17.5 Msi, 
£ , 22 = 1.51 Msi, G’i 2 :=:: C , i 3 =0.78 Msi, C’23 = 0.55 
Msi, and V 12 =0.29. The strength allowables are 
X t = 206.0 ksi, X c =206.0 ksi, YJ=8. 9 ksi, r c =17.8 
ksi, S c 12 = 18.3 Ksi, and S c i3= S C 23=4.8 ksi. Ex- 
perimental results for this test are reported in 
Reference [15]. 

The panel geometry, boundary conditions 
and loading are shown in Figure 5. The panel 
finite element model consists of 6561 nodes and 
6400 elements. Measured geometric imperfection 
from a typical test specimen was included in the 
model. Three integration points through each ply 
thickness are used in the analysis for computation 
of section properties. 

The results for the curved panel loaded in 
compression are shown in Figure 6, where the 
load is plotted against the end shortening dis- 
placement. The filled symbols representing FV1, 
FV2, and FV3 indicate the load level at which 
damage described by field variable FV 1 through 
FV 3 are initiated. The analysis and expermental 
results are in good agreement. The final failure 
load from the experiment is 51.25 kips, which is 
11 % less than the final failure load obtained from 
progressive failure analysis of about 56.8 kips. All 
damage modes initiate after the panel buckling 
and attaining its final load bearing capability. 
Matrix cracking (FVT) initiated at a load level 
of 43.31 kips just after buckling. Fiber-matrix 
shear (FV 2) and fiber failure (FV r 3) initiated at 
a load level of 37.53 kips and a load level of 38.32 
kips, respectively. 

Figure 7(a) through 7(c) show fringe plots 
for matrix cracking (FVT), fiber matrix shear 
(FV2), and fiber failure (FV3) in the bottom ply 
after final failure at an applied displacement of 
approximately 0.1 -in.. Damage initiated within 
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than the experimental failure load. Damage initi- 
ation starts as fiber-matrix shear {FV 2) and fiber 
{FV 3) failure at the same load level of 28.34 kips. 
Matrix cracking {FV 1) is initiated at a load level 
of 30.22 kips. Hence the curved panel can carry 
an additional load of about 11.0 kips (about 27% 
of the final experimental load) after initiation of 
fiber-matrix shear and fiber failure. After matrix 
cracking has been initiated, the panel can carry 
an additional load of 9.0 kips (about 23% of the 
final experimental load). 

All the failure modes initiated near the edge 
of the cutout. Figure 12 shows a fringe plot for 
matrix cracking {FV 1) in the bottom ply after 
the final failure load. These damage locations are 
consistent with experimental observations. How- 
ever experimental observation indicated signifi- 
cant delamination around the cutout. Although 
the initial geometric imperfection are accurately 
represented here, the delamination damage may 
be responsible for the discrepancy between the 
analytical and experimental failure loads with the 
panel exhibiting a catastrophic failure with no 
residual strength. 

Bead-stiffened Panel Loaded in Shear 

The thermoformed bead- stiffened configura- 
tion is an advanced concept for stiffened graphite 
thermoplastic panels. Thermoforming ia a cost 
effective manufacturing method for incorporat- 
ing bead stiffeners. An experimental and analyt- 
ical investigation of these bead-stiffened panels 
was conducted by Rouse [17]. The bead-stiffened 
panels were loaded in in-plane shear loading us- 
ing a picture frame fixture similar to the one de- 
scribed above. It was found that the bead stiff- 
ened panels failed near the curved tip of the stiff- 
ener where large magnitudes of stress resultants 
were predicted. 

The bead-stiffened panel has a laminate 
stacking sequence of [±45/ ± 45/0/ ± 45/90] s 
with a ply thickness of 0.005-in. The mechan- 
ical properties of the material are £u=18.0 
Msi, J E 22 == FfiO Msi, Gi2=(ji 3=C?23 = 0*82 Msi. 
The allowables are X*=300.0 ksi, X c =210.0 ksi, 
Yi=13.0 ksi, y c =31.0 ksi, S c i 2 =27.0 Ksi, and 
Sc 13=S C 23=5.0 ksi. A finite element model of 
the bead-stiffened panel is shown in Figure 13. 
This model consists of 2935 nodes and 2849 ele- 
ments. No geometric imperfection was added to 
the model of the bead-stiffened panel. Three in- 
tegration points through each ply thickness are 
used in the analysis for computation of section 
properties. 


The results for this in-plane shear-loaded 
panel are shown in Figure 14, where the load is 
plotted versus the axial strain (c Xir ) in the top 
and bottom plies (45° ply) at the center of the 
test panel. The curves FV1, FV2, and FV3 in- 
dicate the load level at which damage modes de- 
scribed by the field variable FV1 through F\ 3 
initiate. The analyses results are in good agree- 
ment with the experimental results. The final 
failure load obtained from the experiment is 27.9 
kips, which is 12 % less than the final failure load 
obtained from the progressive failure analysis of 
about 24.4 kips. Damage initiation starts as ma- 
trix cracking {FV 1) at a load level of 10.8 kips, 
whereas fiber-matrix shear {FV 2) and fiber fail- 
ure are initiated at the same load level of 12.8 
kips. Hence the structure can carry an additional 
load of approximately 17.0 kips (60% of the final 
experimental load) after matrix cracking damage 
has been initiated. Even after initiation of fiber- 
matrix shear and fiber failure, the panel can carry 
an additional load of approximately 15.0 kips ( 
53% of the final experimental load). 

Figure 15(a) through 15(c) show fringe plots 
for matrix cracking, fiber-matrix shear, and fiber 
failure in the top ply after the final failure load. 
The damage initiated near the curved tip of the 
bead-stiffener and propagated to the other re- 
gions as shown in these figures. The locations 
for damage are consistent with experimental ob- 
servations. For the case of this stiffened panel, 
initial geometric imperfections are not a critical 
factor in predicting the observed behavior. 

Concluding Remarks 

The results of an analytical and experimen- 
tal study to evaluate the initiation and progres- 
sion of damage in nonlinearly deformed stiffened 
find unstiffened panel are presented. These stud- 
ies are also conducted for panels with cutouts 
and subjected to two loading conditions. The 
progressive failure methodology includes matrix 
cracking, fiber-matrix shear, and fiber failure, 
but ignores delamination failure. The effect of 
initial geometric imperfections is also investi- 
gated as part of the study. 

For a flat panel loaded in in-plane shear load- 
ing, the three failure modes considered in the 
study accurately represent the damage senario in 
the postbuckling regime. The analytically deter- 
mined response, failure modes and damage loca- 
tions compare well with the experimental results 
when the initial geometric imperfection is in- 
cluded in the analysis in a simple manner. When 
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a cutout is introduced, however, delamination oc- 
curs at the hole boundary as an additional fail- 
ure mode resulting in some discrepancies with the 
observed behavior. The analysis results also pre- 
dict a residual strength with the final failure load 
being approximately 10% greater than the exper- 
imental failure load. 

The response of curved panels with and 
without cutouts are studied when loaded in com- 
pression. With measured geometric imperfection 
included in both curved panel models, the re- 
sponse of the panel without cutout compares well 
with experimental results. No delamination oc- 
curs for this case and the failure modes consid- 
ered in this paper develop after buckling. For 
the curved panel with a cutout, however delam- 
ination does occur with the panel failing catas- 
trophically. Unlike the panel without a cutout, 
this panel with a cutout exhibits no residual 
strength. It may be important to include de- 
lamination failure mode to predict the residual 
strength of curved panels with a cutout of the 
type considered here. 

The bead-stiffened panel response, failure 
modes and damage locations are well predicted 
by the analysis results. The residual strength val- 
ues from the analysis and experiment are within 
12 percent with no geometric imperfection in- 
cluded in the analysis. 
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Table 1: Dependence of material elastic 
properties on the field variables 


No 

Matrix 

Fiber-matrix. 

Fiber 

failure 

cracking 

shear 

failure 

£ n 

£n 

£n 

Eu ~>0 

E22 

E22 ~^o 

E22 

E22 — >0 

^12 

V\2 — >0 

0 

t 

1^12 — »0 

G 12 

G 12 

G\2 ~ ^0 

G 12 ->0 

G 13 

G\ 3 

G 13 —>0 

G i3 —>0 

G23 

G23 

G 23 

G23 ~ ►o 

FI' 1=0 

FV'1=1 

FV \= 0 

FV'1=0 

FV'2=0 

FV'2=0 

FV 2=1 

FV'2=0 

FV 2=0 

FV'3=0 

FV'3=0 

FV 3=1 



Figure 1 Schematic diagram of picture frame 
test fixture. 



Figure 2 Finite element model of a flat 
composite panel in the picture frame test 
fixture. 
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Figure 3 Load vs. strain component normal to 
fiber direction in top and bottom plies of flat 
panel loaded in the picture frame test fixture. 



in the top ply of the flat panel. 



Figure 5 Geometry and boundary conditions for 
curved composite panel. 



Figure 6 Load vs. end-shortening displacement 
results for curved panel loaded in compression. 



Figure 7(a) Fringe plot for matrix cracking in 
the bottom ply for curved panel loaded in 
compression. 



Figure 7(b) Fringe plot for fiber-matrix shear in 
the bottom ply for curved panel loaded in 
compression. 
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Figure 7(c) Fringe plot for fiber failure in the 
bottom ply for curved panel loaded in 
compression. 



Figure 8 Finite element model of flat panel with 
a circular cutout. 



Figure 9 Load vs. strain component normal to 
fiber direction in the top and bottom plies of 
flat panel with a cutout loaded in the picture 
frame test fixture. 



Figure 10 Fringe plot for fiber failure in the top 
ply for the flat panel with a cutout loaded in 
the nirtnre frame test fixture. 



Figure 11 Load versus end-shortening 
displacement for a curved panel with cutout 
(d/W= 0.2) and subjected to axial compression 
loading. 



Figure 12: Plot for matrix cracking damage in 
the bottom ply for curved panel with cutout 
(d/W= 0.2) and loaded in compression. 
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Figure 13: Finite element model of the 
bead-stiffened panel in picture frame test 
fixture. 


Load (lbs) 



Figure 14: Load versus axial strain, (e xx ) results 
at the center of bead-stiffened panel. 



Figure 15(b) Fringe plot for fiber-matrix shear 
in the top ply of the bead-stiffened panel 



Figure 15(c) Fringe plot for fiber failure in the 
top ply of the bead-stiffened panel 



Figure 15(a) Fringe plot for matrix cracking in 
the top ply of the bead-stiffened panel 
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